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We propose a "DFT+dispersion" treatment which avoids double counting of dispersion terms by 
deriving the dispersion-free density functional theory (DFT) interaction energy and combining it 
with DFT-based dispersion. The formalism involves self-consistent polarization of DFT monomers 
restrained by the exclusion principle via the Pauli blockade technique. Any exchange-correlation 
potential can be used within monomers, but only the exchange operates between them. The ap- 
plications to rare-gas dimers, ion-rare gas interactions and hydrogen bonds demonstrate that the 
interaction energies agree with benchmark values. 
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The applicability of the density functional the- 
ory (DFT) to calculations of intermolecular potentials 
of van der Waals complexes depends upon a seamless in- 
clusion of the dispersion energy, a long-range correlation 
effect, in the DFT treatment. This goal has been pursued 
vigorously along many lines [1-7] with varying success. 
One promising avenue consists of using an a posteriori 
dispersion correction added to supermolecular DFT cal- 
culations of interaction energy, in the spirit of Ahlrichs 
et al. [8]. For this strategy to succeed two elements are 
necessary: an accurate, nonempirical description of dis- 
persion energy between DFT monomers and a sensible 
dispersion-free description of supermolecular interaction 
energy within DFT. 

A highly accurate and computationally efficient for- 
mulation of dispersion energy is now available from the 
time-dependent DFT as proposed by Misquitta et al. 
[9] and Hefielmann and Jansen [10], referred to as cou- 
pled Kohn-Sham (CKS) dispersion. The second element 
of this strategy, however, has not been available up to 
now. In contrast to the Hartree-Fock (HF) interaction 
energy which is well defined and contains physically in- 
terpretablc effects, the analogous DFT interaction en- 
ergy has neither of these characteristics. Depending on 
a particular functional it may include a variety of ob- 
scure terms. In particular, for all current generalized 
gradient approximation (GGA) formalisms no dispersion 
contribution is accounted for at large intemolecular sep- 
arations R and only some part appears at intermediate 
and small R. As a consequence, approximate exchange- 
correlation (xc) functionals often exhibit an artifactual 
behavior in the long range as well as in the van der Waals 
minimum region. A rigorous DFT+dispersion approach 
should be based on a DFT interaction energy that a pri- 
ori avoids these residual dispersion terms but allows for 
accurate mutual exchange and polarization effects — an 
analog of the HF interaction energy at the DFT level 
of theory. Such a DFT interaction energy could then 



be combined with a dispersion component obtained from 
CKS or other formalisms [11, 12]. 

To this end we adapt the formalism of Pauli-blockade 
Hartree-Fock (PB HF) [13] combined with the Afunc- 
tional subsystem formulation of DFT of Rajchel et al. [14] 
Specifically, the energy of the complex is evaluated from 
the classic Heitler-London formula which takes the an- 
tisymmetrized product of participating monomer wave- 
functions while the monomers are described with the 
KS orbitals from DFT calculations. In the second step, 
one iteratively evaluates the interaction energy between 
two DFT monomers, described by KS determinants, in 
a manner analogous to the HF method. That is, the 
monomers are polarized in each other's fields until self- 
consistency under the constraint of the Pauli exclusion 
principle between monomers. Within the monomers any 
exchange-correlation DFT potential may be employed, 
whereas between monomers the exact exchange potential 
is used to avoid the dispersion contribution. In the third 
step, the dispersion component is a posteriori added. In 
such a way the erratic behavior of approximate xc func- 
tionals is eliminated and all physically important effects 
are included without the problem of the dispersion dou- 
ble counting. Note that the concept of the separated 
monomers resembles the idea of the "range-separation" 
approach in DFT methodology [15]. 

In the present Letter we present an outline of the ap- 
plication of the PB scheme to the calculation of the in- 
teraction energy of DFT monomers. More details on 
derivation and implementation will be published else- 
where. The procedure starts with KS solutions for the 
isolated monomers A and B. The solution is obtained in a 
self-consistent way involving coupled equations for both 
subsystems. The equation for monomer A reads: 
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and its monomer B counterpart is simply generated 
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through the exchange of A and B indices. In Eq. (1), 
[nl superscripts denote iteration numbers, 
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is the standard KS operator (e.g., see Ref. 16) built from 

the monomer A orbitals , the two following terms 

are electrostatic and exchange HF potentials, and the 
last term is the so-called penalty operator enforcing the 
proper mutual orthogonality between A and B monomer 
orbitals so that the intersystem Pauli exclusion principle 
is fulfilled. r\ > is a parameter that does not affect 
the final solution. Monomer A density in terms of KS 
orbitals is 
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The xc potential in (2) requires the asymptotic cor- 
rection of Griming et al. [17] in order to yield 
appropriate interaction energies. The orbitals are 
symmetrically orthogonalized after each iteration, 

yielding ,(&L n+1] j set. The tilde 
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sign (e.g. E[J>) denotes quantities calculated with such 
orbitals. The interaction energy at each iteration is de- 
composed into several terms: 
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~ Ex [Px\ terms are a result of the 
imposition of the intersystem Pauli exclusion principle 
and contain repulsion energy. The expressions for elec- 
trostatic and exchange energies are very simple due to 
the orthogonality of orbitals: 
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Coupling between the subsystems via Eq. (1) leads to the 
system energy lowering, and is referred to as deformation 
energy: 
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represents the final PB energy, calculated with (4) 
using self-consistent orbitals satisfying (1). It results 
from the mutual electric polarization of DFT monomers, 
and, owing to the Pauli blockade procedure, it contains 
exchange contributions. It is related to the induction 
terms of the SAPT formalism with two important ad- 
vantages over the latter: PB gathers all electric polar- 
ization terms to infinity and accounts for accompanying 
exchange effects in a consistent manner within the DFT 
formalism. 

In the original HF-based formulation, the PB proce- 
dure simply restores the supermolecular HF interaction 
energy, and obviously neglects any kind of electron cor- 
relation. For the DFT analog, both monomers are de- 
scribed with the full KS operator, but are coupled us- 
ing HF coulomb and exchange operators [ug lst and v^ ch 
in Eq. (1), respectively] built from KS orbitals. Such an 
approach accounts for intramonomer local electron corre- 
lation, leaving out the intermonomer nonlocal contribu- 
tions. The (fj^ 13 represents then "non-dispersion" part of 
the interaction energy that includes the electrostatic, ex- 
change, and induction components. For rare gas dimers it 
is purely repulsive. In Fig. 1 is compared with the 
supermolecular counterpoise corrected DFT interaction 
energy for Ar2. Calculations employed three DFT func- 




where v\ c is the potential due to monomer A nuclei 
and the intermonomer nuclear-nuclear repulsion 

term. The zero-iteration energy (4) is the DFT ana- 
log of the well-known Heitler-London interaction en- 
ergy: 25™^ = <^n t L (cf. also Cybulski and Seversen 
[18]). The HL energy contains the well defined elec- 
trostatic and exchange interaction contributions, in this 
case between unperturbed isolated DFT monomers. It is 
closely related to the first order energy in the symmetry- 
based perturbation theory based on DFT SAPT (DFT). 



FIG. 1: (Color online) Comparison of PB and supermolecular 
DFT interaction energies for Ax2. 



tionals of the GGA PBE hierarchy, which systematically 
improve on the description of monomers: PBEX [19] (ex- 
change only), PBEREV [20] (local exchange plus corre- 
lation), and PBEO [21] (a hybrid of local and exact ex- 
change plus correlation). All calculations used aug-cc- 
pVTZ basis sets. As seen in Fig. 1 is indeed purely 
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repulsive for Ar2. By contrast, the supermolecular DFT 
interaction energies reveal minima of an unknown origin. 

The total interaction energy, termed PBD for Pauli 
blockade plus dispersion, is obtained by adding to 
the dispersion component obtained from CKS or other 
accurate techniques [11, 12]: 

E i^ U = 4nf + -^disp- (8) 

In this work we use the second-order CKS dispersion 
components combined with the exchange-dispersion term 
from SAPT(DFT) [9, 22, 23]. 

D can now be compared with SAPT(DFT) theory. 
SAPT(DFT) interaction energy 

E™ = jB (i) +E (2) +feF (9) 

includes first-order (electrostatic and exchange) and 
second-order (induction, dispersion and their exchange 
counterparts) contributions derived from KS orbitals. 
<5hf denotes approximate correction for higher-order in- 
duction effects derived at the HF level of theory [24]. 
The SAPT(DFT) terms were calculated using MOLPRO 
program [25]. 

To test the efficiency of the PBD approach, we per- 
formed calculations for several diatomic systems com- 
posed of closed-shell atoms and ions: Ar 2 , ArNa + , and 
ArCl - . Again, the same set of PBE functionals has been 
used with the same aug-cc-pVTZ basis set. The com- 
parison of energies calculated with Eqs. (8) and (9) is 
presented in Figs. 2 through 4. 
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FIG. 2: (Color online) Comparison of PBD and SAPT<5 in- 
teraction energies for At2- 

The results are compared with high-level ab initio 
data from Refs. 26 (Ar 2 ), 27 (ArNa+), and 28 (ArCl - ). 
One can see a remarkably good agreement of PBD and 
benchmark data, similar to that of SAPT5, for Ar 2 , and 
ArNa + . The agreement for ArCl - is somewhat worse, 
however for this system the benchmark calculation may 
be up to 10 % too shallow. 




FIG. 3: (Color online) Comparison of PBD and SAPT<5 in- 
teraction energies for ArNa + . 




-4 1 ' ' ' ' ' ' 1 

5.5 6 6.5 7 7.5 8 8.5 9 

FIG. 4: (Color online) Comparison of PBD and SAPT<5 in- 
teraction energies for ArCl - . 



Another test of the proposed approach was performed 
for hydrogen-bonded and other molecular complexes. 
The results for a representative set of these complexes 
are shown in Table I. The PBD interaction energies 
are compared with SAPT^ and with benchmark values. 
PBD and SAPT<5 share the same values of dispersion 
energy. The benchmarks correspond to basis set satu- 
rated CCSD(T) results. The comparison indicates that 
PBD leads to very reasonable interaction energies for 
hydrogen-bonded systems. As mentioned above, the two 
components <£^ t L and #™ have a clear physical interpre- 
tation of electrostatic-plus-cxchange interaction of two 
unperturbed monomers and mutual monomer polariza- 
tion contribution restrained by exchange, respectively. It 
should also be noted that whereas SAPTi5 provides re- 
sults of equally high quality, it is dependent upon inclu- 
sion of (5hf which is substantial. 
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TABLE I: Comparison of interaction energies and their com- 
ponents from PBD calculations with SAPT(DFT) and bench- 
mark values (in mH). Equilibrium geometries are from indi- 
cated references. PBEO functionals with aug-cc-pVTZ basis 
set are used for monomers (aug-cc-pVDZ for CO2). 



System S^~ a E^ Dc S^ 3 E^ 1 " 1 '* Benchmark 

^812 -9.00 [29] 

-6.48 -7.22 [30] 

-3.16 -3.10 [30] 

-4.78 -5.05 [29] 

-13.2 -13.55 [30] 

-0.948 -0.813 [29] 

-2.11 -1.80 e [31] 



(H 2 0) 2 -1.02 -4.94 -8.56 -1.45 

(HF) 2 -0.693 -4.78 -7.59 -1.19 

(HG1) 2 1.65 -0.548 -3.48 -0.976 

(NH 3 ) 2 -0.479 -2.03 -5.34 -0.522 

H2O-HF -1.15 -10.2 -14.8 -3.02 

(CH 4 ) 2 0.705 0.631 -0.852 -0.0378 

(CQ 2 ) 2 0.85 0.501 -1.74 -0.0673 



"Electrostatic + exchange. 
'Electrostatic + exchange 
c Electrostatic + exchange 
''Higher-order induction. 

e CCSD(T) result for the global minimum from [31]. 



mutual polarization, 
mutual polarization 



dispersion. 



In summary, we presented a new treatment of interac- 
tion energy between two DFT monomers which is exactly 
dispersion-free. It contains two physically-meaningful 
terms: HL interaction energy and self-consistent polar- 
ization which is restrained by the exact exchange. This 
interaction energy which is an analog of the supermolecu- 
lar HF interaction energy can be combined with existing, 
reliable treatments of dispersion to yield the first theoret- 
ically sound DFT+dispersion approach. The computa- 
tional cost of our approach is essentially the same as that 
of the routine KS calculations, and it scales with the size 
and number of fragments as the standard KS equations. 
It is worthwhile to note that for the test cases of this Let- 
ter the results proved to be remarkably insensitive to the 
asymptotic correction — an issue which deserves further 
studies. 
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